Chain metallicity and antiferro-paramagnetism competition in underdoped 
YBa 2 Cu 3 06+ s : a first principles description 
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We describe from advanced first principles calculations the energetics of oxygen doping and its 
relation to insulator-metal transitions in underdoped YBa2Cu306+:r. We find a strong tendency 
of doping oxygens to order into non-magnetic Cu 1+ Oa; chains at any x. Ordering produces one- 
dimensional metallic bands, while configurations with non-aligned oxygens are insulating. The 
Cu 2+ C>2 planes remain insulating and antiferromagnetic up to a threshold between £=0.25 and 
0.5, above which a paramagnetic normal-metal state prevails. The in-plane antiferro-paramagnetic 
competition depends on x, but only weakly on the ordering state of the chains. 

PACS numbers: Valid Pacs 



A central puzzle in the physics of the archetypal high- 
T c superconducting cuprate YE^CusOg+a; is the oc- 
currence of several seemingly concurrent phase transi- 
tions -structural, magnetic, order-disorder, insulating- 
metal, and superconducting- in the underdoped region 
of the (T,a;) phase diagram. Superconductivity sets on 
at £ c ^0.35, so that an insulator-metal transition (IMT) 
should occur around the same doping. This IMT, proba- 
bly associated with a change in magnetic state, is thought 
be related to an order-disorder transition in the Cu-0 
chains holes form a Fermi glass in the low-doping 
Mott insulator up to a density h~0.07 per CuC>2 (roughly 
corresponding to x c ), above which they start percolating 
in the planes leading to a conventional metal state. Ex- 
periments, however, suggest that CuO chains undergo 
an independent IMT at lower doping than the CuC>2 
planes. At strong to moderate underdoping (a;~0.2-0.8), 
infrared experiments identify the electromagnetic re- 
sponse of the chains as that of a Tomonaga-Luttinger 
liquid-like one-dimensional conductor, although the chain 
fragments are too short (^15-400 A) to support dc cur- 
rents across macroscopic regions (except for weak un- 
derdoping (£=0.95) whereby they are conducting [!, H|). 
Optical data [H| show that oxygens cluster in metallic is- 
lands with local orthorhombic symmetry even at minute 
doping, and that metallic percolation through the chains 
occurs at a doping as low as £=0.1-0.2. 

It is fair to say that the relation of chain ordering 
with one-dimensional metallization and with the IMT 
in the Cu02 planes at higher doping has never been 
described thoroughly from first principles. Here, us- 
ing advanced first principles electronic structure calcu- 
lations, we a) study the connection between intra-chain 
one-dimensional metallicity and oxygen ordering in (non- 
magnetic) CuO chains as function of doping, and b) an- 
alyze the competition of the antiferromagnetic (AF) in- 
sulating and paramagnetic (PM) metallic states of the 
Cu02 planes in different chain ordering configurations. 
The picture we arrive at is as follows. YBa2Cii306 is, 



as expected, a Mott-Hubbard insulator. At small x it 
remains everywhere insulating inasmuch as doping oxy- 
gens are disorderly distributed through the chains. A 
dispersed, one-dimensional metallic band appears when 
oxygens orderly align into CuO chains. Chain alignment 
is found to be always favored, consistently with the evi- 
dence that chain nucleation starts immediately at x>0 
and that metallic percolation in the chains is already 
detectable at x^O.l 0, @|. The dispersed band has a 
strongly dominant chain-atom character, and holes re- 
main largely bound to the doping oxygen 0(1) and the 
nearby Cu(l) ions. Only a minor hole fraction is trans- 
ferred to the Cu02 planes (indicating a rather sharp de- 
coupling of the two playgrounds): the ensuing small vari- 
ations of orbital occupations govern the very tight com- 
petition between the in-plane insulating AF state and 
the conventional metallic PM phase. The AF state gives 
way to the PM metallic state at a doping between x=0.25 
and 0.5, consistent with the overall metallization thresh- 
old implied by the superconductive onset: intriguingly, 
this crossover occurs independently of the way 0(1) 's are 
distributed in the chains, i.e. on the order-disorder inter- 
play governing the in-chain IMT. 

The realistic description of underdoped cuprates is 
troublesome for first-principles band theories such as 
local-spin density (LSDA) or generalized gradient ap- 
proximation (GGA) density-functional theory Q, which 
incorrectly describe the antiferromagnetic Mott insulat- 
ing state as non-magnetic and metallic, and therefore 
are not usable to study the IMT. Here we circumvent 
this issue employing the pseudo-self-interaction correc- 
tion method (PSIC) [7], which isparticularly accurate 
for strongly correlated materials [8J. As a backup and 
countercheck, we also use GGA+U Q and GGA calcu- 
lations. PSIC calculations are carried out with a plane 
wave basis ultrasoft pseudopotential [l(| method, with 
30 Ryd cut-off energies. For GGA and GGA+U we em- 



ploy the PAW method as implemented in the VASP [11 1 



code. We use 6x6x6 special k-point grids for total energy 
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FIG. 1: OR-DOS of AF YBa 2 Cu 3 06. 0(2) and 0(4) are 
in-plane and apical oxygens, respectively. Only the most im- 
portant orbitals (d for Cu, Ba, and Y, and p for O) are shown; 
solid line is the sum of the 3 nearly degenerate t 2fl DOS. For 
Op p y ~Pz, while for Ot Px~p y - Supercell: 48-atoms 2x2. 
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FIG. 2: OR-DOS of AF YBa 2 Cu 3 06.25 for the broken-chain 
configuration (a) in FigO Labels are same as FigQ] 



calculations, 250 special k-points and linear tetrahedron 
interpolation method for DOS. The GGA+U parameters 
are U=10 eV, J=l eV, giving the best match with PSIC 
bands. At different doping concentrations, we use su- 
percells ranging from 24 ( V2 x \/2x 1 primitive cells) up 
to 96 atoms (4x2x1). Unless otherwise specified, the 
results shown below are obtained by PSIC. 

In FigQ] we show the orbital-resolved density of states 
(OR-DOS) for tetragonal AF YBa 2 Cu 3 6 , which shows 
a Mott-Hubbard charge-transfer band gap ~1.2 eV (the 
observed photoconductive threshold at ~1.5 eV jjj). The 
valence band top (VBT) and conduction band bottom 
(CBB) are hybrids of unpolarized O (p^, p y ) and, re- 
spectively, majority (occupied) and minority (unoccu- 
pied) Cu d x 2_ y 2 states. The fundamental transition only 
involves orbitals within Cu02 planes; transitions involv- 
ing final states with apical 0(4) p z and chain Cu(l) d 2 2 
states (as well as Ba and Y states) start only above ~3.5 
eV (more details will be presented elsewhere). 

Now we consider oxygen inclusion in the chains. In 
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FIG. 3: (Color online) Doping configurations in the Cu(l)-0 
chains. Dark (blue) and light (yellow) circles are O(l) and 
Cu(l) ions, respectively, (a): x = 0.25, cell 2x2; (b) and (c): 
x = 0.25, cell 4x2; (d), (e), and (f): x = 0.5, cell 2x2. 
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FIG. 4: Band structure of AF YBa2Cu3 06.25- Left: broken- 
chain (configuration (a) of Fig[3]). Right: full chains (config- 
uration (b), same Figure). 



FigE] we show the OR-DOS at x=0.25 calculated with 
one doping oxygen per 2x2 unit cell (configuration (a) of 
FigEl) corresponding to an alternated array of parallel, 
isolated Cu(l)-0(1)-Cu(l) segments: this is our broken- 
chain configuration. The main feature induced by 0(1) 
is a spin-degenerate hole state ~1 eV above the VBT, 
strongly localized on 0(1) itself with dominant p^ char- 
acter (x is the Cu(l)-0(1) direction). This hole is resid- 
ually extended (FigJ2) on Cu(l) d x 2_ y 2 (through dom- 
inant pdtr coupling) and, more conspicuously, on Cu(l) 
d z 2 . This is because d z 2 is the highest (most antibonding) 
Cu(l) d state (see Cu(l) OR-DOS in FigHJ). The hole 
further propagates to apical 0(4) p 2 states (see FigEJ , 
and up into the Cu02 planes. The calculated orbital 
decomposition of the hole band gives 32.3% on 0(1), 
33.5% on the two adjacent Cu(l), 27% on the apical 0(4) 
closer to the latter, 2.5% on 0(2), 1.3% on Cu(2), and fi- 
nally 2.3% on Ba. The corresponding band structure (left 
panel of FigH]) shows the flat hole band running through 
the whole Brillouin zone just below the undoped conduc- 
tion band (actually almost touching it at T). 

To evaluate the effects of oxygen alignment at £=0.25 
(configuration (b) of Fig [3]) we employed a 4x2 (i.e. 98- 
atom) supercell, whose band structure is reported in the 
right panel of FigHJ Now the hole band is dispersed and 
crosses the Fermi energy, so the system is metallic. The 
dispersion is one-dimensional: no metallicity is injected 
into the Cu0 2 layer, as signaled by the band flatness 
along the M-L k 2 segment. Even at £=0.5 the OR-DOS 
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FIG. 5: Orbital-resolved DOS of ordered-chain AF 
YBa2Cu306.5- The calculation is done in configuration [d] 
of Fig[3] Note that the while dispersed band crossing Ef 
has sizable Cu(l)-0(l)-0(4) character, the local DOS in the 
Cu02 planes is still insulating. 



of the chain-ordered AF configuration (Figj5]) shows that 
the hole band does not affect the band gap in the AF 
Cu02 layer, and spreads instead on the same in-chain 
orbitals previously seen to be involved in the localized 
hole decomposition. 

So far we have shown that the order-disorder and one- 
dimensional IMT in the chains are directly related. This 
relation holds for £=0.25 and £=0.5, so it must hold 
for any 0.25<x<0.5. Of course the modeling of disorder 
in our calculations is limited by the size of our super- 
cell: for a realistic sample of disordered distributions we 
may envisage a hole energy distribution spanning (and 
eventually closing) the ~1 eV energy interval above the 
VBT. So, rather than as an optical insulator, the chains 
electronic ground state is probably better described as a 
hopping-conductive Fermi glass, in agreement with the 
mainstream interpretation. However, the limits in simu- 
lating disorder effects do not change the conclusion that 

a) only chain- alignment produces hole band dispersion, 

b) the hole band is one-dimensional in character, and the 
corollary that c) as long as the AF ordering is retained, 
metallic conductivity does is not transferred outside the 
chains. 

We now discuss the energetics of chain- aligned and un- 
aligned configurations at £=0.25 and £=0.5. The data 
for both AF and PM ordering, obtained with different 
approaches, are listed in Table HI The chain-ordered con- 
figuration is the most stable by a safe margin, irrespective 
of Cu02 magnetic ordering. Thus AE can be thought as 
the cost of breaking a single chain moving an 0(1) to the 
nearby isolated segment. These results imply a strong 
drive towards chain formation and hence one-dimensional 
metallic behavior. Although disorder does not allow the 
formation of macroscopic chains at low doping 0, our 
calculated energies confirm the observation that chain 
nucleation starts at very low doping, also accompanied 
by local one-dimensional metallic percolation. 

Next, in Figl6] we compare the AF and PM total en- 



TABLE I: Total energies for the O(l) configurations in Fig[3] 
([b] and [c] for x=0.25, [d], [e], and [f] for 1=0.5), and AF and 
PM phases (values labeled with * and ** refers to GGA+U 
and GGA calculations, respectively). 
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ergies for chain-ordered ([b] or [d]) and chain-alternated 
([c] or [f]) configurations as a function of x. The energy 
zero is the lowest phase at each doping. The AF phase 
is favored up to £=0.25 whereas the PM phase is favored 
at and above £=0.5. 
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FIG. 6: Comparison of the energies of AF and PM phases as 
function of x for broken-chain (top) and ordered-chain (bot- 
tom) configurations. 

Three points are worth noticing. First, the AF^PM 
crossover (with attendant in-plane insulating^metallic 
transition) occurs in the doping interval where super- 
conductivity sets on. Second, the AF-PM competition 
is barely affected by the 0(1) distribution within the 
chains; this implies in both instances a weak chain-plane 
coupling with minimal hole transfer to the planes, con- 
firming the high inefficiency of oxygen doping. Third, the 
AF-PM competion is quite tight (a few tenths of meV per 
formula), in line with the tendency of the AF phase to 
form cluster spin-glass states and the coexistence of mi- 
croscopic AF and PM domains [l2| ■ 

The AF-PM competition can be interpreted analyzing 
the changes in orbital occupation upon doping, reported 
in Table ITT1 Each 0(1) takes over ^0.4 electrons. This 
quantity is quite independent of magnetic ordering, but is 
visibly larger for configuration (d) than for (f): full-chain 
order enhances charge transfer between the doping O and 
its surrounding. The vast majority of the hole released by 
0(1) remains in the apical units (Cu(l)-0(l)-0(4)), and 
only a fraction is transferred to the planes. The most 
receptive environment is PM<j, which collects 13.9% of 
the hole charge into the planes, followed by PMf (11.2%), 
AF d (8.7%) and AF f (6.4%). 

We can understand this classification in terms of mag- 
netic and chain ordering. The PM phases have in-plane 
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TABLE II: Differences (in 10 -3 unit charge) in atomic orbital 
occupation between x=0 and £=0.5 for PM and AF order- 
ing and chain configurations (d) and (f) (see Fig(3|. Positive 
(negative) charges are hole (electron) fractions. The number 
of equivalent atoms in the cell is in square brackets. Subscript 
"1" indicates atoms next-neighbors of 0(1) or in apical units 
containing 0(1); subscript "2" non-neighbors of, or in units 
without, 0(1): e.g. 0(4)i is first-neighbor to a Cu(l)-0(1) 
unit; 0(2)i is on-top of 0(1), 0(2)2 and 0(2)3 are co-planar 
to 0(2)i; Cu(l)i and Cu(l)2 are first- or second-neighbors of 
0(1); Cu(2)i and Cu(2) 2 are on-top of Cu(l)i and Cu(l) 2 , 
respectively. 
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CuC>2 bands at Ejr, allowing hole injection in the highly 
receptive Cu(2) atoms. In AF phases the Cu(2) states 
are swept off the E f region and are marginally affected 
by hole injection. Table [IT] shows that the PM in-plane 
hole population is mainly on Cu(2), and in the AF phase 
in- holes mainly resides on 0(2), with Cu(2) having small 
electron-like excess charge. As for 0(1) configuration, (d) 
ordering is a better environments for hole transfer than 
(f). The metallic band originating from chain ordering, 
with its hybridization with apical atoms, is better suited 
to charge transfer than insulating disordered states. 

Finally, the charge transfer efficiency is the key to un- 
derstand the PM-AF competition. Charge injection in 
the planes tends to stabilize the PM phase at increasing 



doping. At £=0.5, the AF-to-PM transition is accompa- 
nied by increase of ~0.05 holes per Cu02 unit for both 
the (d) or (f) configurations, i.e. the AF-PM hole charge 
imbalance does not depend on the 0(1) configuration. 
This explain the substantial similarity of the two panels 
in Fig [51 

In summary, we described from first principles the ef- 
fects of oxygen doping in underdoped YP^CusOe+z- We 
found chain alignment energetically favored with respect 
to broken chains at any doping concentration. Since oxy- 
gen alignment ordering produces one-dimensional metal- 
lic bands, chain metallicity is formally present at any 
doping, although over short distances only. We also con- 
sidered the energy competition of the insulating AF and 
metallic PM state in the Cu02 planes, showing that the 
latter phase is favored above a threshold doping located 
between x=0.25 and £=0.5. The in-plane AF-PM transi- 
tion is caused by the (sparing) chain-to-planc hole trans- 
fer, but is quite independent of the oxygen ordering in 
the chains. 
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